Revisiting the chain magnetism in SruCi^Ozii: Experimental and numerical results 
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We study the magnetism of the hole doped CuC>2 spin chains in Sri4Cu24C>4i by measuring the 
Electron Spin Resonance (ESR) and the static magnetization M in applied magnetic fields up to 
14 T. In this compound, the dimerized ground state and the charge order in the chains are well 
established. Our experimental data suggest that at low temperatures the Curie-like increase of 
M as well as the occurrence of the related ESR signal are due to a small amount of paramagnetic 
centers which are not extrinsic defects but rather unpaired Cu spins in the chain. These observations 
qualitatively confirm recent ab initio calculations of the ground state properties of the CuC>2 chains 
in Sri4Cu24 04i. Our complementary quantum statistical simulations yield that the temperature 
and field dependence of the magnetization can be well described by an effective Heisenberg model 
in which the ground state configuration is composed of spin dimers, trimers, and monomers. 
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I. INTRODUCTION 

Hole doped antiferromagnets with an inhomogeneous 
ground state have attracted much attention in recent 
years. A prominent example is provided by high tem- 
perature (high-T c ) superconductors in which low dimen- 
sionality, s = 1/2 quantum magnetism, and charge de- 
grees of freedom are intimately related. In particular, a 
nonuniform distribution of charges in the Cu-0 planes of 
the high-T c layered cuprates seems to play an important 
role in establishing the superconductivity^ The quasi 
one-dimensional (ID) edge-sharing Cu02 chain com- 
pounds such as (Sr, Ca, La)i4Cu2404i, Nai +x Cu02, or 
Ca2+ x Y2- x Cu50io constitute another class of cuprates 
with charge inhomogeneitiesiSi 3 ^^ The interplay of spin 
and charge degrees of freedom in one dimension, as is 
evidenced by both theoretical and experimental results, 
is even more pronounced as compared to the layered 2D 
situation, yielding unusual ground states and novel exci- 
tations both in the spin and charge sectors. 

In this work we focus on the low-temperature mag- 
netic properties of the so-called "telephone number" com- 
pound Sri4Cu2404i, where a composite structure with 
CU2O3 spin ladders and CUO2 spin chains is realized. 
Magnetically, the ladders do not contribute significantly 
at low temperature due to a large spin gap of A ~ 
3gQj^ 6j7j8j9 a remarkable property of this compound is 
self-doping, i.e. the material is intrinsically hole doped 
with 25% of holes per Cu site. As most of the holes are lo- 
calized in the chains their concentration in this subunit is 
very high, close to 60% per Cu in the chain formula unit 



(c.f.u.). This influences strongly the magnetism of the 
spin chain resulting in an inhomogeneous ground state. 
Two Cu spins which are separated by a hole site cou- 
ple antiferromagnetically (AF). If they are surrounded 
by two localized holes at both sides they form almost in- 
dependent dimers. Such a dimerized nonmagnetic state 
in the chain subsystem of Sri4Cu2404i is experimentally 
well establishedii2iiiii2ii&ii The occurrence of charge or- 
der, which is an important ingredient of the dimer sce- 
nario in Sri4Cu240 4 i, is supported e.g. by NMR and 
ESR dataji5ii& In a very recent numerical studjiii Gelle 
and Lepetit (GL) highlight that structural modulations 
due to the pseudoperiodicity of the chains and the lad- 
ders might be crucial for the electron localization and for 
the occurrence of the dimer state. Remarkably, these cal- 
culations suggest that the ground state is not perfectly 
dimerized even in the limiting case where all holes are lo- 
cated in the chains. Ideally, in that case, where six out of 
the ten sites in the chain unit cell are holes, the remain- 
ing four spins could be coupled to two AF dimers with 
no free spins left. In contrast to that, GL arrive at the 
surprising result that it is energetically favorable if a few 
spins are left unpaired. Then the stable configuration 
corresponds to 1.7 dimers and 0.5 free spins per c.f.u., 
respectively. Therefore, the occurrence of a Curie contri- 
bution to the magnetic susceptibility of Sri 4 Cu2404i at 
low temperatures, which has originally been attributed to 
the free spins, might be an intrinsic property of the hole 
doped Heisenberg Cu-0 spin chain due to the distortion 
by structural modulations. 

The aim of our paper is thus twofold. First, we 
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show experimental magnetization and electron spin res- 
onance (ESR) data which give strong indications that 
a small amount of free spins, that gives rise to a fi- 
nite magnetization and to an ESR signal at low tem- 
peratures, is not extrinsic but mainly reside in the chain 
subunit of Sri4Cu2404i. This finding qualitatively sup- 
ports the model calculations by GL in Ref. Sec- 
ondly, we present our own quantum statistical simula- 
tions which rely on the numerical diagonalization of an 
effective Heisenberg Hamiltonianii2iiS This Hamiltonian, 
which depends parametrically on hole positions and on 
the screened Coulomb interaction between them, also 
yields an inhomogeneous ground state configuration of 
spins and holes. It turns out that the ground state is 
dominantly built of the aforementioned dimers, but also 
contains weakly coupled monomers and/or trimers which 
not only leads to the Curie-like behavior of the magne- 
tization at low temperatures, but also yields a nearly 
perfect quantitative description of the low-temperature 
magnetic response of the Cu spin chain. The very good 
agreement also suggests that other arrangements of spins 
and holes are well gaped from the ground state configu- 
ration. 
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II. EXPERIMENTAL RESULTS 

In this section we present ESR and magnetiza- 
tion measurements on a high quality single crystal of 
Sri4Cu2404i. For ESR experiments an X-Band Bruker 
ESR spectrometer has been used, which operates at a fre- 
quency of 9.5 GHz. The magnetization data have been 
collected with a homemade vibrating sample magnetome- 
ter (VSM) in external magnetic fields up to 14 T and 
with a SQUID magnetometer in a field of 1 T. The sin- 
gle crystal grown by the floating zone technique^ has 
been previously thoroughly characterized by measuring 
the magnetization, the thermal expansion^ the electri- 
cal and thermal transport the inelastic neutron scat- 
tering cross section^, and the ESR response^ 

In order to recall the magnetic properties of 
Sr 14 Cu240 4 i we show in Fig. [J a) the temperature de- 
pendence of the magnetization M(T) of our single crys- 
tal in a magnetic field of B = 1 T applied parallel to 
the crystallographic c-axis (chain direction). The M(T) 
dependence can be very well fitted with Eq. Q compris- 
ing two terms, the magnetization Mdimer of AF coupled 
dimers in a concentration Nd with an exchange constant 
J/k B = —134 K, and the Curie magnetization Mcurie 
owing to a small concentration Ng of free spins (see, e.g. 
Ref. ED 
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The very good agreement of the experimental M (T) 
curve with Eq. give evidence of the dimerized ground 



FIG. 1: (color online) Magnetization of Sri4Cu2404i versus 
temperature for a constant magnetic field of B = 1 T (a) and 
versus applied magnetic field for 2.5 K< T < 20 K (b). B 
was applied along the c-axis. Lines are fits to the data (see 
text). 



state which in addition has been proven for our sample 
by inelastic neutron scattering and thermal expansion 
measurementsii2i2i In Fig. IHb) we present the magneti- 
zation versus applied magnetic field, M(B), at tempera- 
tures of 20 K, 15 K, 10 K, 4.2 K, and 2.5 K, respectively. 
At all temperatures the dependence of M on B can be 
nicely fitted as a sum of the temperature independent 
contribution xo'B owing to diamagnetism and Van-Vleck 
paramagnetism of the Cu ions and a Brillouin function 
B s= i describing the magnetization of free s = 1/2 spins 
with a concentration N$ 



M(B) = xo 



B. 



gii B {B + \M) 
2k B T 



(2) 

With a g factor of g c — 2.04 (see below) one obtains 
Xo ~ 1 • 10~ 5 emu/Mol Cu, the mean field parameter 
A = and Ns — 0.01/Cu. We note that the magnetiza- 
tion at T = 2.5 K saturates in a field of B — 14 T at a 
very small value of about 1 • 10 -2 fi B /Cu justifying that 
the magnetic response at low temperature is caused by 
about 1% of free s — 1/2 spins while there is no response 
from the remaining 99% of the spins. Moreover, we em- 
phasize that the determination of N$ from our high field 
magnetization data is much more straightforward and 
accurate as compared to the conventional analysis of the 
low temperature contribution to M(T) using the Mcurie 
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term in Eq. (JTJ. Therefore, precise knowledge of N$ from 
the analysis of M(B) curves allows for an accurate de- 
termination of the number of dimers No- Substituting 
Ns = 0.01/Cu in the Mcurie term in Eq. one ob- 
tains from the fit to the experimental M(T) dependence 
No = 0.0738/Cu. Thus the number of magnetic Cu sites 
amounts to Ns + 2Nd ks 3.78/f.u. in our Sri4Cu2404i 
single crystal which differs from a previous less precise 
estimate that employed the analysis of M(T) curves of 
polycrystalline samples^ 

We note that the concentrations of dimers and free 
spins, Njj and N$, respectively, are significant observ- 
ables in the recent ab initio calculations of the low energy 
electronic properties of the Cu-0 chain in Sri4Cu2404i 
by GLiii Therefore, the accurate determination of Nfj 
and Ns from our magnetization data provides an impor- 
tant check for the relevance of the GL model with respect 
to the low energy physics of Sri4Cu2404i. Comparison of 
our data with the numerical results shows that the num- 
ber of dimers Njy found in our analysis roughly agrees 
with the calculations of GL for the case that all holes are 
localized in the chains, N§ L = 0.0708/Cu . The number 
of free spins Ng, however, is smaller than theoretically 
predicted N^ L = 0.021/Cu by a factor of two. This is 
surprising because the experimental value Ns comprises 
not only the response owing to the imperfect dimeriza- 
tion of the spin chain, but also should include the re- 
sponse of paramagnetic defects which may occur in real 
crystals. Thus, theoretical calculations should yield a 
smaller number of free spins than experimental data. Be- 
sides the quantitative discrepancy this comparison sug- 
gests that, at least in the framework of the GL model, 
the free spins in a concentration Ns — 0.01/Cu found 
in our Sri4Cu2404i single crystal are not extrinsic para- 
magnetic defects but mainly reside in the spin chain^ 
This suggestion is corroborated with the analysis of the 
anisotropy of the ESR signal and of the magnetization. 

Detailed studies of the ESR response from the spin 
chains in the family of the "telephon e num ber" com- 
pounds have been published in Refs. Il6l25i Here we 
present the ESR data of the Sri4Cu2404i single crystal 
in the low temperature regime. Below 20 K the ESR sig- 
nal is mainly due to the free spins, which also contribute 
to the Curie-like increase of the magnetization, compare 
Fig-H^a). At higher temperatures the ESR response is 
dominated by the dimer spins which are thermally acti- 
vated. The evolution of the g factor of the ESR signal in 
the crossover regime around 20 K is shown in Fig. EI The 
g factor of the Cu spins in a dimer is anisotropic with the 
values g c = 2.045 for the magnetic field B parallel to the 
chain direction ( c-axis) and gt = 2.284 for the magnetic 
field B perpendicular to the plane of the chains (&-axis), 
respectively. These g factors are temperature indepen- 
dent. The anisotropy of the g tensor is determined by 
the crystal field splitting of the Cu orbital state of e g 
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FIG. 2: Temperature dependence of the g-factor along c- and 
6-axis, respectively, as obtained from ESR data. Dashed lines 
extrapolate g of the dimers. Below ~20 K, the response is due 
to small deviations from the perfectly dimerized state. 



and tig symmetry 

8 • A „ 2- A 

9c =2-- ,g b = 2-— . (3) 

Here A = —0.1 eV is the spin-orbit coupling constant 
and A yZtZX and A xy are the energy differences between 
relevant orbital states^ The anisotropy of the g tensor 
above ~ 20 K is temperature independent and is typical 
for a Cu ion in an approximately square planar coordi- 
nation of the oxygen ligandsji§*2i i.e. the coordination 
occurring in the chain of Sri4Cu2404i. Remarkably, be- 
low 20 K the g factor for both orientations changes only 
slightly giving strong indication that the ESR response at 
low Ts comes from a paramagnetic Cu ion having a very 
similar environment to that of the chain Cu ions. A small 
change in the anisotropy of the g factors can be explained 
by a slightly different charge distribution around a free 
Cu spin which has in the chain two neighboring holes, 
both to the left and to the right hand side, whereas a Cu 
spin in a dimer has only one neighboring hole on the one 
side and two holes on the other side, compare also Fig. [3] 

The fact that the free spins exhibit the similar 
anisotropy as the dimerized ones is confirmed by our 
magnetization data. Our analysis of M(T) shows that 
the ratio of the Curie constants which describe the re- 
sponse of the free spins along b- and c-axis amounts to 
Cb/C c = g 2 /g 2 — 1.25. This ratio agrees well with our 
ESR data. This fact is illustrated by Fig. [3] in which 
M/g 2 is shown for B parallel to the 6-axis and the c-axis, 
respectively. After being corrected for the anisotropy of 
the g factor the magnetization data for two orientations 
coincide almost perfectly in the whole temperature range. 
Note that in our analysis we have also considered the 
anisotropy of the Van-Vleck magnetism which is related 
to the anisotropy of the g factors asSi 
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FIG. 3: (color online) Magnetisation divided by the p-factor 
and corrected by anisotropic Van-Vleck magnetism along the 
respective axes of Sri4Cu24 04i for a constant magnetic field 
B — 1 T versus temperature. 
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To summarize these results, our experimental data 
strongly suggest that the free spins which cause the 
Curie-like response in Sri4Cu2404i are not extrinsic 
paramagnetic defects but reside in a not perfectly dimer- 
ized chain. Thus our experimental study qualitatively 
confirms the results of GL in Ref. \TA which predict the 
occurrence of such free spins in the hole doped chain 
as a result of structural modulations. As noted above, 
the theoretical results suggest, however, a larger Curie- 
like response than observed experimentally. In addition, 
while GL assume slightly more than 60% of holes in the 
chains of Sri4Cu2404i, recent NEXAFS experiments im- 
ply that there are less holes in the chains^ i.e. there is a 
larger number of magnetically active spins in the chains. 

To resolve this contradiction we perform a numerical 
study to be presented in the next section which gives ev- 
idence for trimers in spin chains with less than 60% hole 
doping. Trimers consist of three next-nearest neighbor 
spins with a hole between each spin, which are confined 
by two holes at each side, compare Fig. |SJ The pres- 
ence of trimers might lift the discrepancy between the 
calculated number of paramagnetic defects in the dimer- 
ized state and the smaller experimental value of Ns- For 
small external fields, trimers respond similarly to free 
spins. The number of magnetically active spins involved 
in this response is, however, three times larger. The first 
non-trivial excitation of trimers occurs at 1.5 J ~ 200 K, 
which can not be observed in our M(B) data at 2.5 K be- 
cause this is far above the accessible field range. On the 
other hand, our M(T) data do not rule out this scenario 
because at higher temperatures there are indeed addi- 
tional contributions to M which are usually ascribed to 



the magnetism of the ladders and to the melting of the 
charge order. Our numerical analysis in section lTTll shows . 
that trimers (instead of or in addition to single spins) 
perfectly explain our experimental data. In contrast, the 
presence of quadrumers, which exhibit the spin gap of 
0.659 J ~ 90 K is excluded by our M(T) data shown in 
Fig-H^a). Our data also rule out that the Curie- like re- 
sponse is due to nearest neighbor dimers, trimers, etc., 
which are coupled ferromagnetically. 



III. EXACT DIAGONALIZATION STUDIES 

The dynamics of a hole-doped chain system consti- 
tutes a real challenge for theoretical investigations. Espe- 
cially the evaluation of thermodynamic quantities both 
as function of temperature and magnetic field is pro- 
hibitively complicated even for moderate system sizes. In 
the following we are using an effective spin Hamiltonian 
which depends parametrically on hole positions^ This 
ansatz is similar to a simple Born-Oppenheimer descrip- 
tion where the electronic Hamiltonian (here spin Hamil- 
tonian) depends parametrically on the positions of the 
classical nuclei (here hole positions). Each configura- 
tion c of holes and spins defines a Hilbert space which 
is orthogonal to all Hilbert spaces arising from different 
configurations. The Hamilton operator H(cj of a certain 
configuration c is of Heisenberg type, i. e. 
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Juv (c) are the respective exchange parameters which de- 
pend on the configuration of holes. J < describes an- 
tiferromagnetic coupling, J > ferromagnetic coupling. 
For the theoretical results presented in this article three 
exchange parameters are used, see Fig. The strongest 
and antiferromagnetic exchange J = —134 K is across 
one hole. The exchange accross two holes Jy = 15 K is 
ferromagnetic as is the exchange Jnn — 100 K of neigh- 
boring spins. Periodic boundary conditions are applied 
for the following calculations, i.e. large rings are consid- 
ered instead of chains with open boundaries. 

It is obvious that different configurations of spins and 
holes should also be energetically different. Besides their 
different magnetic ground state energies resulting from 
they also differ by their Coulomb interaction be- 
tween the holes and possibly by the Coulomb interac- 
tion with the local environment. In the following we 
take the electrostatic hole-hole repulsion by means of a 
screened Coulomb potential into account. r n = 2.75 A is 
the distance between nearest neighbor sites on the ring. 
The accurate value of the dielectric constant e r is un- 
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FIG. 4: Exchange parameters used for the exact diagonaliza- 
tion: J = —134 K, J« = 15 K, and Jjvjv = 100 K. 



known. Several attempts have been undertaken to es- 
timate the dielectric constant which yielded values for 
e r up to 30. 27 - 29 - 30 If the interaction with the local - 
modulatedii^i^SiS - environment is neglected, it is nec- 
essary to assume a rather small dielectric constant e r J$ 3 
in order to simulate the magnetization dataM> Under 
such conditions only the energetically lowest-lying spin- 
hole configuration c contributes to the magnetization for 
T < 200 K. 
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FIG. 5: Building blocks of the hole doped Cu02 chains: open 
circles denote hole sites, arrows denote spin sites. The first 
block contains a single spin (s), the second a dimer (d), and 
the third block hosts three spins (t), which form a short chain. 
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FIG. 6: (color online) Magnetization versus temperature for 
B — IT: Experimental data (B || c-axis) are given by crosses. 
The results of a complete numerical diagonalization are de- 
picted by a solid curve for iV s = 17 and Nh — 26 as well as 
by a dashed curve for N s = 17 and Nh = 25. 



be investigated. To this end Hamiltonian JBJ is numer- 
ically completely diagonalized for the currently largest 
possible odd number of spins N s = 17. The number 
of holes is chosen as JVh = 25 or N-^ — 26 in order to 
yield about 60 % holes on the chain. Periodic bound- 
ary conditions are applied. It turns out that the ground 
state configurations are sequences of the building blocks 
shown in Fig. [3] In the case of N s = 17 spins and 
N& = 26 holes the resulting ground state sequence is 
Ci7,26 = sdddddddd, i.e. it contains only dimers except 
for one single spin block. This ground state configuration 
is very similar to the one found using density functional 
theory calculationsii2i2£ In the case of N s = 17 spins and 
— 25 holes the chain consists of dimers and one block 
of three spins, i.e. it has the sequence 017,25 = tddddddd. 
In both cases other configurations of spins and holes are 
energetically well separated. 

Having determined all energy eigenvalues and mag- 
netic quantum numbers, the magnetization can be eval- 
uated for the two ground state configurations. Figure 
shows the resulting magnetization curves as a function 
of temperature for an applied field of B = 1 T. The 
solid curve displays the magnetization curve for N s = 17 
and A^h = 26, i.e. configuration cfj.726 = sdddddddd, 
whereas the dashed curve shows the result for N s = 17 
and Nix = 25, i.e. 017,25 = tddddddd. Both curves are 
rather close to the experimental data given by crosses. 

The magnetization has also been determined as a func- 
tion of applied field for two temperatures T = 2.5 K and 
T = 4.2 K. Figure [7| shows that both configurations, i.e. 
Ci7,26 = sdddddddd and 017,25 = tddddddd reproduce the 
data with high accuracy. 
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FIG. 7: (color online) Magnetization versus applied field for 
T = 2.5 K and T = 4.2 K: Experimental data (B || c-axis) 
are given by symbols. The results of a complete numerical 
diagonalization are depicted by a solid curves for N s = 17 
and Nh = 26 as well as by dashed curves for N„ = 17 and 
iV h = 25. 



IV. CONCLUSION 



The assumption of a perfectly dimerized chain was 
studied in Ref. In the following imperfect chains will 



In this paper we present experimental as well as numer- 
ical evidence that the Curie-like contribution to the low- 
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temperature magnetization of Sri4Cu2404i is a genuine 
property of the CuC>2 chain subsystem. The hole-doped 
chains are not regular sequences of weekly interacting 
spin dimers, but instead they host a small percentage 
of almost free spins and/or spin trimers. This finding is 
in accord with recent density functional calculationspiZ*^ 
Nevertheless, our own simulational studies favor the oc- 
currence of trimers instead of single spins which would 
be consitent with an intrinsic hole doping of the chains 
of less than 60 % as was inferred from recent NEXAFS 



experiments 
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